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Abstract 

We propose a 2-dimensional cellular automaton model to simulate pedestrian traffic. 
It is a Vmax = 1 model with exclusion statistics and parallel dynamics. Long-range 
interactions between the pedestrians are mediated by a so called floor field which 
modifies the transition rates to neighbouring cells. This field, which can be discrete 
or continuous, is subject to difi'usion and decay. Furthermore it can be modified by 
the motion of the pedestrians. Therefore the model uses an idea similar to chemo- 
taxis, but with pedestrians following a virtual rather than a chemical trace. Our 
main goal is to show that the introduction of such a fioor field is sufficient to model 
collective eff'ects and self-organization encountered in pedestrian dynamics, e.g. lane 
formation in counterfiow through a large corridor. As an application we also present 
simulations of the evacuation of a large room with reduced visibility, e.g. due to fail- 
ure of lights or smoke. 
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1 Introduction 



Considerable research has been done on the topic of traffic flow using methods 
from physics during the last decade [1-8]. Cellular automata inspired by the 
pioneering works [9,11,10] compose by now an important class of models. Most 
studies have been devoted to one-dimensional systems, where several analytic 
approaches exist to calculate or approximate the stationary state. 
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The majority of these models deals with particles which can move by more 
than one cell per time step (maximal velocity v^ax > !)• Furthermore, it seems 
to be widely accepted that the most suitable update procedure is the parallel 
(synchronous) update. Both open and periodic boundary conditions have been 
considered, where problems with open boundaries are generally harder to treat 
analytically (for a review, see [12]). 

On the other hand, pedestrian dynamics has not been studied as extensively 
as vehicular traffic, especially using a cellular automata approach. One reason 
is probably its generically two-dimensional nature. In recent years, continuum 
models have been most successful in modelling pedestrian dynamics. An im- 
portant example are the social force models (see e.g. [4,8,13] and references 
therein). Here pedestrians are treated as particles subject to long-ranged[^ 
forces induced by the social behaviour of the individuals. This leads to (cou- 
pled) equations of motion similar to Newtonian mechanics. There are, however, 
important differences since e.g. in general the third law ("actio = reactio") is 
not fulfilled. 

In contrast to the social force models our approach is closer in spirit to the 
general strategy of modelling (elementary) forces on a microscopic level by 
the exchange of mediating particles which are bosons. It is therefore similar 
to active walker models [14,15] used so far mainly to describe trail formation, 
chemotaxis (see [16] for a review) etc. Here the walker leaves a trace by modi- 
fying the underground on his path. This modification is real in the sense that 
it could be measured in principle. For trail formation vegetation is destroyed 
by the walker and in chemotaxis he leaves a chemical trace. In contrast in our 
model the trace is virtual. Its main purpose is to transform effects of long- 
ranged interactions (e.g. following people walking some distance ahead) into a 
local interaction (with the "trace"). This allows e.g. for a much more efficient 
simulation on a computer. 

Cellular automata for pedestrian dynamics have been proposed in [17-22]. 
These models can be considered as generalizations of the Biham-Middleton- 
Levine model for city traffic [10]. Most works have focussed on the occurrence 
of a jamming transition as the density of pedestrians is increased. All mod- 
els have Vmax = 1, except for the generalization proposed in [22] which is 
used for analyzing evacuation processes on-board passenger ships. The other 
models use a kind of "sublattice-dynamics" which distinguishes between differ- 
ent types of pedestrians according to their preferred walking direction. Such 
an update is not easy to generalize to more complex situations where the 
walking direction can change. To our knowledge so far most other collective 
effects encountered empirically [4,23,24,26,25] have not been reproduced using 
these models. Another discrete model has been proposed earlier by Gipps and 
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Marksjos [27]. This model is somewhat closer in spirit to our model than the 
cellular automata approaches of [17-22] since the transitions are determined 
by the occupancies of the neighbouring cells. However, also this model can not 
reproduce all the collective effects. In [28] a discretized version of the social 
force model has been introduced. The repulsive potentials by the pedestrians 
are stored in a global potential, with pedestrians reacting to the gradients 
of this global potential. Although this model is able to reproduce collective 
effects it suffers from some drawbacks [29]. It is not flexible enough to treat 
individual reactions to other pedestrians, and collision-avoidance is not always 
guaranteed for velocities greater than 1. 

First we discuss some general principles we took into account in the devel- 
opment of our model. In contrast to vehicular traffic the time needed for 
acceleration and braking is negligible. The velocity distribution of pedestrians 
is sharply peaked [23]. These facts naturally lead to a model with Vmax = 1 
(if space is discrete), i.e. only transitions to nearest neighbours^ are allowed. 
Furthermore, a greater Vmax would be harder to implement in 2 dimensions, 
especially when combined with parallel dynamics, and reduce the computa- 
tional efficiency. The number of possible target cells increases quadratically 
with the interaction range. Furthermore one has to check whether the path 
is blocked by other pedestrians. This might even be ambigious for diagonal 
motion and crossing trajectories. 

To keep the model simple, we strongly emphasize the principle to provide the 
particles with as little intelligence as possible and to achieve the formation of 
complex structures and collective effects by means of self-organization. Effec- 
tively, there is absolutely no intelligence (i.e. look-ahead distances or multiple 
moves per update step depending on the distribution of occupied neighbour 
sites) in our model. In contrast to older approaches we do not make detailed 
assumptions about the human behaviour. Nevertheless the model is able to 
reproduce many of the basic phenomena. 

The key feature to substitute individual intelligence is the floor field. Apart 
from the occupation number each cell carries an additional quantity (field) 
which can be either discrete or continuous. This field can have its own dy- 
namics given by diffusion and decay coefficients. 

Interactions between pedestrians are repulsive for short distances. One likes 
to keep a minimal distance to others in order to avoid bumping into them. In 
the simplest version of our model this is taken into account through hard-core 
repulsion which prevents multiple occupation of the cells. For longer distances 
the interaction is often attractive. E.g. when walking in a crowded area it is 
usually advantageous to follow directly behind the predecessor. Large crowds 
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may also be attractive due to curiosity. 

In order to produce a flow around obstacles in a simple way we also present 
a variant of our model where the pedestrians can be in one of two modes (or 
moods), "happy" or "unhappy". These two modes are distinguished by their 
fluctuations. " Happy" pedestrians try to move in a preferred direction whereas 
"unhappy" pedestrians move in a more random fashion. This is sufficient to 
avoid a jamming transition due to obstacles at unrealistical low densities. 

With two particle species moving in opposite directions, each with its own 
floor field, effects can be observed which are so far only achieved by continuous 
models: lane formation and oscillation of the direction of flow at doors. We 
consider this model to be another proof of the ability of cellular automata to 
create complex behaviour out of simple rules and the great applicability to all 
kinds of traffic flow problems. 

The model can be used together with models for route selection which assign 

certain routes to each pedestrian. It only assumes that at every timestep for 
each pedestrian a transition matrix (matrix of preferences) is given. 



2 Model 



The underlying structure is a 2-dimensional grid which can be closed period- 
ically in one or both directions. Each cell can either be empty or occupied 
by exactly one particle (pedestrian). The size of a cell corresponds to approx- 
imately 40 X 40 cm^. This is the typical space occupied by a pedestrian in 
a dense crowd [26]. For special situations it might be desirable to use a finer 
discretization, e.g. such that each pedestrian occupies four cells instead of one. 
In this paper, however, we concentrate on the simplest case which seems to be 
sufficient for most purposes. The update is done in parallel for all particles. 
This introduces a timescale into the dynamics which can roughly be identi- 
fied with the reaction time treac- In the deterministic limit, corresponding to 
the maximal possible walking velocity in our model, a single pedestrian (not 
interacting with others) moves with a velocity of one cell per timestep, i.e. 
40 cm per timestep. Empirically the average velocity of a pedestrian is about 
1.3 m/s [26]. This gives an estimate for the real time corresponding to one 
timestep in our model. It is approximately 0.3 sec, i.e. of the order of the 
reaction time treac, and thus is consistent with our microscopic rules. 
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Fig. 1. A particle, its possible transitions and the associated matrix of preference 

M={M,j). 



2.1 Basic Rules 



Each particle is given a direction of preference. ^Prom this direction, a 3 x 3 
matrix of preferences is constructed which contains the probabilities for a move 
of the particle. The central element describes the probability for the particle 
not to move at all, the remaining 8 correspond to a move to the neighbouring 
cells. The probabilities can be related to the velocity and the longitudinal 
and transversal standard deviations (see appendix A for details). In practice 
all particles of the same species share the values of these parameters and in 
consequence the same matrix. In the simplest case the pedestrian is allowed 
to move in one direction only without fluctuations and in the corresponding 
matrix of preference only one element is one and all others are zero (see Fig. 1). 



This ansatz can easily be extended by fixing the direction of preference for each 
cell separately, e.g. to handle structures inside buildings. Then the particles 
would use the matrix belonging to the cell they occupy at a given step. 

In each update step, for each particle a desired move is chosen according to 
these probabilities. This is done in parallel for all particles. If the target cell is 
occupied, the particle does not move. If it is not occupied and no other particle 
targets the same cell, the move is executed. If more than one particle share 
the same target cell, one is chosen according to the relative probabilities with 
which each particle chose their target. This particle moves while its rivals for 
the same target keep their position (see Fig. 2). 

The matrix used to achieve maximum flow is clearly the simplest case of uni- 
directional fluctuation-free motion described above. Transversal fluctuations 
reduce the flow by introducing interference between lanes. However, this set- 
ting is not sufficient to avoid e.g. jams behind obstacles. To escape a jammed 
situation, the particles need a mode where they can select between different 
cells to move backwards or sideways in order to eventually make their way 
around the obstacle. 
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Fig. 2. Solving conflicts according to the relative probabilities for the case of two 
particles with matrices of preference M^^^ and M^^^ . 

The rules presented up to here are a straightforward generalization of the CA 
rules used so far for the description of traffic flow [17-21]. The main difference 
is that in principle transitions in all directions are possible and each pedestrian 
j might have her own preferred direction of motion characterized by a matrix 
of preferences M^^^ . The only interaction between particles taken into account 
so far is hard-core exclusion. 



2.2 Floor Field 



In order to reproduce certain collective phenomena it is necessary to introduce 
further longer-ranged interactions. In some continuous models this is done 
using the idea of a social force [4,8,13]. Here we introduce a different approach. 
Since we want to keep the model as simple as possible we try to avoid using a 
long-range interaction explicitly. Instead we introduce the concept of a floor 
field which is modified by the pedestrians and which in turn modifies the 
transition probabilities. This allows to take into account interactions between 
pedestrians and the geometry of the system (building) in a unified and simple 
way without loosing the advantages of local transition rules. The floor fleld 
modifies the transition probabilities in such a way that a motion into the 
direction of larger fields is preferred. 

The floor fleld can be thought of as a second grid of cells underlying the 
grid of cells occupied by the pedestrians. It can be discrete or continuous. In 
this paper we will give examples for both variants. In general we distinguish 
between static and dynamic floor flelds. The static floor field S does not evolve 
with time and is not changed by the presence of pedestrians. Such a fleld can be 
used to specify regions of space which are more attractive, e.g. an emergency 
exit (see the example in Sec. 3.1) or shop windows. This has an effect similar 
to a position-dependent matrix of preference but is much easier to realize. 
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In contrast the dynamic floor field D is modified by the presence of pedestrians 
and has its own dynamics, i.e. diffusion and decay. Usually the dynamic floor 
field is used to model a ("long-ranged") attractive interaction between the 
pedestrians. Each pedestrian leaves a "trace", i.e. the floor field of occupied 
cells is increased. Since the total transition probability is proportional to the 
dynamic floor field it becomes more attractive to follow in the footsteps of 
other pedestrians. Explicit examples where such an interaction is relevant will 
be given in Sec. 3. The dynamic floor field is also subject to diffusion and 
decay which leads to a dilution and finally the vanishing of the trace after 
some time. 

In general the transition probability pij in direction {i,j) (see Fig. 1) is given 

byQ 

Pij = NMij Dij Sij{l-nij). ( 1 ) 

Here Uij is the occupation number of the target cell in direction (2, j), i.e. Uij = 
for an empty cell and Uij = 1 for an occupied cell. Therefore transitons to 
occupied cells are forbidden. is a normalization factor to ensure J2{i,j) Pij = 1 
where the sum is over the nine possible target cells. In Sec. 2.4.2 we will also use 
a slightly different form for the transition probabilities which is more general 
than (1). 

The update rules of the full model including the interaction with the floor 
fields then have the following structure: 

(1) The dynamic floor field D is modified according to its diffusion and decay 
rules (see Sec. 2.3 and 2.4). 

(2) For each pedestrian, the transition probabilities for a move to an unoccu- 
pied neighbour cell {i,j) is determined by the matrix of preferences and 
the local dynamic and static floor fields, e.g. pij oc MijDijSij. 

(3) Each pedestrian chooses a target cell based on the probabilities of the 
transition matrix P = {pij). 

(4) The conflicts arising by any two or more pedestrians attempting to move 
to the same target cell are resolved, e.g. using the procedure described in 
Sec. 2.1. 

(5) The pedestrians which are allowed to move execute their step. 

(6) The pedestrians alter the dynamic floor field of the cell they occupied 
before the move. 

The explicit form of the interaction between the pedestrians and the floor field 
and the dynamics of the floor field will be specified in Sec. 2.3 and 2.4. 



Note that this is not a product of matrices but just the product of the corre- 
sponding matrix elements. 
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If more than one pedestrian species exists (i.e. two groups moving in opposite 
directions), each species interacts with its own floor field. In the simplest case 
these fields are independent from each other. 



2.3 Discrete floor fields 

In the discrete case the fields are realized through noninteracting particles 
which do not obey a hard-core exclusion principle. Therefore they will be called 
bosons in the following. Since the particles corresponding to the pedestrians 

are not allowed to share a cell these will be called fermions. The fermions 
couple to the bosons locally which drives the fermions in a preferred direction 
and induces a long-range interaction between the fermions. 

The essence of our approach is that for each fermion the probability to jump 
into a direction with a larger number of bosons is increased. Thus, the motion 
is simply driven by gradients in the floor field, i.e. in the density of the bosons. 

The first type of bosons (s-bosons) is completely static. At the beginning of 
a simulation for each cell {x^y) the occupation number of s-bosons Ts{x,y) 
is fixed to a specific value. Furthermore, at the beginning every cell is void 
of bosons of the second type, the dynamic bosons (d-bosons). Whenever a 
fermion jumps from site {x,y) to one of the neighbouring cells, the d-boson 
occupation number of cell {x,y) is increased by one (fermions leave a trace): 

rd{x,y) ^Td{x,y) + 1. (2) 

After all motions of the fermions during one timestep have been performed, 
the oldest d-boson of each cell is destroyed with probabihty a, if the lifetime 
of this boson is larger than one (i.e. it has been created during the previous 
update step or earlier). 

Now that the two bosonic floor flelds have been introduced, the update proce- 
dure for the fermions can be given. At every discrete time step t ^ t + 1 each 
fermion verifies which of its neighbouring cells are empty {uij = 0). The 
transition probability to occupied neighbouring cells is set to zero. Thus, the 
probability for a jump from the center cell (0, 0) to an unoccupied neighbour 
site is given by 

Pij = Nexp{pJsAs{iJ)) ■ exp(/3JdAd(i, j)) • (1 - riij) ■ dij (3) 

where 

A,{i,j)^r,ii,j)-rsiO,0) and Aa{i, j) ^ r^ii, j) - rdiO,0). (4) 
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N is again a normalization factor to ensure J2{ij)Pij = 1- di is a correction 
factor taking into account the direction the particle in cell has been com- 
ing from. The variables Jj and control the coupling strength between the 
fermions and the s-bosons and the d-bosons, respectively. P plays the role of an 
inverse temperature. Note that the d-bosons lead to a long range interaction 
between fermions in space and time. 

The correction factor is introduced in order to prevent that the fermions are 
not confused by their own trace. One has to distinguish between three cases: 
If the fermion at has been sitting at i at time t — 1, the d-boson sitting on 
top of i has been left by the fermion under consideration. Setting 

dij = ex.p{-(3Jd) (5) 

this fermion is not taken into account in the calculation of the transition 
probabilities. If during the time step from t — 1 ^ t the fermion has moved 
into the direction of the vector pointing from (0, 0) — > this motion shall 

be enhanced 

= exp(/3Jo). (6) 

Therefore, Jq is a parameter which can be used to tune the inertia of the 
fermions. In all other case d^j equals one. 

This prescription is not free of collisions. Therefore, if m fermions try to per- 
form a move onto the same site, only one of these fermions is allowed to 
perform this move. This fermion is picked at random with probability l/m. 
Of course one can also use the method described in Sec. 2.1 for the resolution 
of the conflicts. For the problem studied in Sec. 3.1 the details of the conflict 
resolution turned out to play no important role and we therefore used the 
simpler rule. 

2.4 Continuous floor field 

In the continuous variant each cell j of the floor field carries a continuous 
field value fj between and 1. The basic purpose of the floor field is again to 
determine the transition probabilities of the pedestrians. 

In the example studied in Sec. 3.2 we will only use a dynamic floor field, but a 
generalization which includes also a static field is straightforward. Since we are 

interested in applications related to the flow around obstacles we introduce 
two kinds of states in which pedestrian can be: "happy" or "unhappy". A 
pedestrian becomes " unhappy" if several consecutive desired moves could not 
be carried out due to conflicts. She then changes her strategy, i.e. her matrix 
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of preferences (see Sec. 2.4.2). The interaction of the pedestrians with the floor 
field is then as follows: "Happy" pedestrians locally increase the field, and a 
large field aids "unhappy" pedestrians to become "happy" again, in a sense to 
be specified below. This is sufficient to produce a flow around obstacles, e.g. 
lane formation (see Sec 3.2) or oscillations of the direction of flow at doors. 

Without the distinction between the two states pedestrians would have the 
tendency to pile up in front of obstacles. If a constant flow from behind exists 
it will become increasingly difficult for these pedestrians to turn around and 
avoid the obstacle and therefore the pile will grow. For static obstacles this 
effect is related to the fact that our pedestrians have minimal intelligence. 
This assumption should be well justified in situations where they move in 
unknown territory and at a reduced visibility, e.g. due to smoke or failing 
lights. In normal situations the pedestrians can see a static obstacle (e.g. a 
wall) from some distance and will try to avoid them early. This can easily 
be incorporated in our model using a static fioor field which becomes smaller 
just in front of the obstacle and thus reduces the corresponding transition 
probabilities. For dynamical obstacles, e.g. other pedestrians moving in the 
opposite direction, one way of avoiding unrealistic jamming properties is the 
introduction of different modes. 



2.4-1 Diffusion and decay 

The dynamic fioor field F is subject to diffusion and decay. It evolves according 
to 

OF 

— = D-AF-6-F (7) 

dt ^ ^ 

which is discretized in the standard manner. Here D is the diffusion constant 
and S the decay constant. The ranges are restricted to D e [0, |] and S e [0, |] 
to suppress oscillations in the fioor fields and insure that the values do not 
leave the interval [0, 1]. 



2.4-2 Floor field affects the pedestrian's desire and state 

As mentioned above there are different ways to model the interaction of the 
fioor fields with the pedestrians. Here we present a preliminary solution for 
the situation of pedestrians with minimal intelligence which can probably be 
simplified further. The form of the transition rates is slightly more general than 
(1). Furthermore we allow the pedestrians to be in two different modes (or 
moods) described by different matrices of preference. The transitions between 
these two modes are controlled by the ffoor field. The rates are, however, 
subject to two restrictions: 
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• A uniform floor field should not alter the matrix of preferences. 

• A non-uniform floor field should be able to change a matrix element from 
zero to a nonzero value. 

A pcr-clement addition of the matrix of preferences and the floor field violates 
the first principle, while a multiplication by e.g. Boltzmann factors violates 
the second. Therefore, we are working with a compromise between the two by 
slightly generalizing eq. (1) 

Pij = N {Mij + 62) • exp ( {Fij - Favg) ■ h) , (8) 

where F denotes the floor fleld matrix, F^vg the floor fleld averaged over all nine 
relevant cells and P = (pij) the transition matrix which has to be normalized 
by a normalization factor N. The parameters of this rule are 61 and 62- 

We introduce a second mode by switching to a different matrix of preferences. 
The mode described so far will be called happy, while the unhappy mode is 
realized by a different matrix of preferences which is simply characterized by 
greater standard deviations and a reduced velocity. This happens typically in 
high-density situations. The motion then becomes less directed and the larger 
fluctuations help to avoid clogging. Happy pedestrians which could not move 
to their desired target fleld in several consecutive timesteps (the exact number 
is a parameter of the model, in our simulations set to 3) enter the unhappy 
mode, while unhappy pedestrians become happy again after a certain number 
of consecutive allowed desires (4 in our simulations). The optimal values for 
these two parameters depend on the preferences of the given simulation (i.e. 
maximum flow versus flexible obstacle avoidance). Of course this requires a 
minimal per-pedestrian memory consisting of two counters. 

One could easily alter this deflnition by allowing a continuous choice of matrix 
by interpolating between the two states. This would correspond to introducing 
a continuous spectrum of moods instead of the discrete states "happy" and 
"unhappy". 

In addition to the mechanism for transitions between the modes described 
above, a unhappy pedestrian changes to a happy one immediately if the value 
of the floor field at its cell of origin is greater than a certain threshold. This 
leads to a smooth integration into a calm region of high fiow of pedestrians 
which have run into an obstacle or a jam. 

2.4-3 Pedestrians affect the floor field 

After a pedestrian has completed a certain number of total allowed moves 
(the value of which is usually 3) in the happy mode, the value of the floor fleld 
of its originating cell is increased. This introduces a third counter residing in 



11 



the per-pedestrian memory. The counter is set to zero at mode changes. The 
prescription to alter the floor fleld reads 

F^F + min((l-F)^i, ^2). (9) 

The parameters of this rule arc gi G [0, 1] and g2 G [0, 1]. The way how pedes- 
trians affect the floor fleld can certainly be altered slightly without changing 
the overall behaviour. 

It is important to ignore this change in the field while modifying the matrix in 
the next update step to avoid artifacts (pedestrians moving backwards without 
reason) . 

With this extension, the jam behind a single obstacle composed by a line 
of several forbidden cells can be reduced significantly. The trade-off lies in an 
overall reduced ffow, as the unhappy pedestrians gravely disturb the previously 
unhindered movement of the pedestrians which pass to the sides of the obsta- 
cle. This evokes the need to find out for each pedestrian individually whether 
it might be necessary to switch to the happy state instantaneously, depending 
on the local situation. This should be achieved without the introduction of 
per-pedestrian intelligence. 



3 Simulations 



In the following we describe the results of simulations of two typical situations, 
i.e. the evacuation of a large room [25] (e.g. in the case of a flre) and the 
formation of lanes in a large corridor [13]. We use different variants of the 
basic model in order to elucidate the potential of the different approaches. 



3. 1 Evacuation of a large room 



For simphcity, in the case of discrete floor flelds pedestrians are only allowed 
to move in north [N), west {W), south [S), and east {E) direction, which 
leads to the following form of the matrix of preferences: 



M 



^ Q Mjv ^ 

Mw Mo Me 
\ Ms ) 



(10) 



This choice means no severe restriction since transitions into the diagonal 
directions can be implemented quite easily. 
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In our simulations we have investigated the behaviour of people leaving a 
quadratic room with one door only. The s-bosonic field has been chosen such 
that the occupation number of s-bosons decreases radially from a maximum 
value at the door to zero at the corners opposite to the door. Typical stages 




Fig. 3. People leaving a room with one door only. Displayed are three typical stages 
of the dynamics. 

of the dynamics are shown in Fig. 3. 



As an example we have studied the infiuence of the lifetime of (i-bosons (i.e. 
their decay probability a) on the evacuation time, i.e., the time it takes for all 
people to leave the room. We have seen that if the coupling strength Jg to the 
static bosons is rather large, the evacuation time increases with an increase of 
the lifetime of (i-bosons (see Fig. 4). Most interestingly, if becomes smaller, 
the best evacuation times are found when the lifetime of rf-bosons is fixed at 
some intermediate values (see Fig. 5). 
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Fig. 4. Mean evacuation times Tg^ as function of the decay probability a of the 
d-bosons. The parameters are: = 2, = Jq = 1, /? = 10. The errorbars display 
the mean standard deviation. 



This finding is very interesting, because it has the consequence that the at- 
tractive interaction of particles can lead to crucial differences in the particles' 
behaviour. These changes become more severe if the particles have no clear 



13 



2200 
2000 
1800 
1600 
1400 
1200 
1000 
800 
600 
400 
200 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 O.S 



Fig. 5. Same as Fig. 4, but with Jg = 1/2. All other parameters are the same. 

idea what the best way to the next exit is. In addition, one can see that fluc- 
tuations become much more dominant if Jg goes to zero. Therefore, in case of 
evacuation simulations, studying average evacuation times or - even worse - 
looking at one sample only might lead to wrong conclusions. 



3.2 Lane formation in a long corridor 



We present simulations of a rectangular corridor which is populated by two 
species of pedestrians moving in opposite directions. Parallel to the direction of 
motion we assume the existence of walls. Orthogonal to the direction of motion 
we investigated both periodic and open boundary conditions. The length of 
the corridor is set to 200 cells. Widths of 15 and 25 cells have been used. 

With periodic boundary conditions, the density of pedestrians is fixed for each 
run. The program ensures that the overall number of pedestrians is evenly 
divided by the numbers for the different species. For open boundaries, wie fix 
the rate at which pedestrians enter the system at the boundaries (ASEP style 
insertion rates). The pedestrians leave the system as soon as they reach the 
opposite end of the corridor. 

This model clearly provides the option for a complete jam. The jamming prob- 
ability with periodic boundaries at constant density increases with the length 
of the system. An open system can be thought of as the limit of an infinitely 
long periodic system, although density and entry rates do not correspond 
absolutely (the density in the open system is always higher than twice the 
insertion rate). 

The update rules have the same structure as described in Sec. 2.2. Only step 
(6) is modified to 



14 



(6') The pedestrians change their mode if necessary based on their history 
and the floor flcld. They alter the (dynamic) floor fleld of the cell they 
occupied before the move. 



We performed several runs for different densities and insertion rates, respec- 
tively. The focus of our attention is the parameter range where the transition 
from a stable flow to a complete jam takes place. The complete set of param- 
eters for the simulations can be found in appendix B. 



□ Cellular Automaton 
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Fig. 6. Snapshot of a simulation with p = 0.17, w = h = 50. The left part shows 
the parameter control. The central window is the corridor and the light and dark 
squares are right- and left-moving pedestrians, respectively. The right part shows 
the floor fields for the two species. 

Figure 6 shows the graphical frontend running a simulation of a small periodic 
system. The lanes can be spotted easily, both in the main window showing 
the cell contents and the small windows on the right showing the floor fleld 
intensity for the two species. 
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Fig. 7. Velocity profile of a periodic system with p = 0.10. 
To obtain information about the lanes we accumulated the pedestrian veloc- 
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Fig. 8. Velocity profile of an open system taken at x = ^ with a = 0.04. 
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Fig. 9. Velocity profile of an open system taken at x = ^ with a = 0.04. 

ities at a cross section perpendicular to the direction of flow. This is done 
according to the formula jn+i = jn ■ f + v, where j is the accumulated value 
and V = 0, 1 is the velocity of the pedestrian crossing the hne. r < 1 is set to 
such a value that the characteristic number of contributing pedestrians is 100. 
Selected profiles are shown in Figs. 7-9. The values of the other parameters 
are given in the table in Appendix B. Qualitatively our results are in good 
agreement with those of [30] where the lane formation has been interpreted as 
an optimal self-organization process. 

It is obvious that the lane formation in the periodic system works far better 
than in the open system. The floor field leads to an effective attraction of 
identical pedestrians while different pedestrian species separate. This results 
in the formation of a stable pattern in the periodic case. 

In a certain density regime, these lanes are metastable. Spontaneous fiuctua- 
tions can disrupt the fiow in one lane causing the pedestrians to spread and 
interfere with other lanes. Eventually the system can run into a jam by this 
mechanism. The average time after which the system is blocked by a jam is 
an interesting observable which depends on the density of pedestrians. We 
observe large fiuctuations of this observable which require many samples to 
find statistically significant information. 
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Wc have also found the formation of an odd number of lanes under certain 
conditions. This corresponds to a spontaneous breaking of the left-right sym- 
metry of the system. 

Due to the complexity of our model, the computational speed is significantly 
lower compared to the original models of traffic flow. A typical value measured 
on a SUN Sparc-10 workstation is 0.24 mega updates per second. It should 
still be way faster than continuous models. 



4 Conclusions 

We have introduced a stochastic cellular automaton to simulate pedestrian 
behaviour. We focused on the general concept and the effects which can be 
observed with the basic approach, i.e. particle attraction and repulsion be- 
tween identical and different particles respectively and lane formation. 

The key mechanism is the introduction of the floor field which acts as a substi- 
tute for pedestrian intelligence and leads to collective phenomena. This floor 
field makes it possible to translate spatial long-ranged interactions into non- 
local interactions in time. The latter can be implemented much more efficiently 
on a computer. Another advantage is an easier treatment of complex geome- 
tries. In models with long-range interactions, e.g. the social-force models, one 
always has to check explicitly whether pedestrians are separated by walls in 
which case there should be no interaction between them. Furthermore, the 
computational effort in these models increases proportionally to the square 
of the number of individuals. In contrast, in our approach it increases only 
linearly with the system size which is usually fixed. 

The general idea in our model is similar to chemotaxis. However, the pedestri- 
ans leave a virtual trace rather than a chemical one. This virtual trace has its 
own dynamics (diffusion and decay) which e.g. restricts the interaction range 
(in time). It is realized through a dynamcical floor fleld which allows to give 
the pedestrians only minimal intelligence and to use local interactions. To- 
gether with the static floor field it offers the possibifity to take different effects 
into account in a unified way, e.g. the social forces between the pedestrians or 
the geometry of the building. 

We presented a rather general form of the model. Not all the features are 
needed in all the cases. E.g. for lane formation we do not need a static field. A 
static field might lead to a "pinning" of lanes. We have shown that our model is 
a good starting point for realistic applications since it is able to reproduce the 
basic phenomena encountered empirically. In contrast the other CA models 
so far have not been shown to exhibit some of the collective phenomena. 
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e.g. lane formation etc. Other features, e.g. oscillations at doors [13], have 
also been observed in our simulations [31]. Quantitative results will presented 
elsewhere. The model can also be applied to more complex geometries and 
various characteristics of a crowd can be simulated without major changes. 
So it should be possible to study the effects of panic (see [25] and references 
therein) . 

The description of pedestrians using a cellular automaton approach allows for 
very high simulation speeds. Therefore, we have the possibility to extract the 
complete statistical properties of our model using Monte Carlo simulations. 
This knowldege is of major importance if one wants to estabhsh risk manage- 
ment techniques that are nowadays used for the hedging of financial assets all 
over the world [32]. 
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A Construction of the matrix of preferences 



The aim of this appendix is to show that the matrix of preferences can be 
directly related to observable quantities, namely the average velocities and 
their fluctuations. The procedure explained here to construct the matrix of 
preferences is not essential for the model. One could freely choose the nine 
matrix elements to achieve the desired behaviour of the pedestrians. However, 
it is not straightforward to choose five independent probabilities (one is de- 
termined by normalization and three by symmetry) in a consistent way. It is 
therefore convenient to look for a simpler principle, which might even simplify 
some calculations. 

We consider first a one-dimensional setup of three adjacent fields which rep- 
resents the velocities in / := { — 1,0,1} from left to right. To these cells the 
probabilities po and pi are assigned. The values of the average velocity v 
and of the standard deviation a are the parameters of this construction. This 
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leads to three conditions which uniquely determine the probabilities: 

Ep^ = 1 (A.l) 

(A.2) 

(A.3) 



Not all combinations of f e [—1,1] and a e [0, 1] are allowed. One finds 

P-i = ^{(^^ + v^- v) (A.4) 
Po = 1 - ((t' + v"") (A.5) 
Pi^l{'J^ + v^ + v) , (A.6) 

where a is confined to the interval [u;, ah] with 



(Ti. 



v\ — 



(A.7) 
(A.8) 



These restrictions are shown in figure A.l. 




Fig. A.l. ai and CTh depending on v. 

To create a 3 x 3 matrix M, one creates two such sets: {p-i-,Pf)-,Pi} with the 
parameters {w, c7„} which correspond to the forward and backward movement 
and {g'-i, Q'o, ?i} with the parameters {0,(Jt} for the (symmetric) transversal 
movement. These are simply multiplied: 



1i ■ Pj 



(A.9) 



which produces a movement to the right for positive v. This matrix is normal- 
ized by construction. We have now achieved a reduction of free parameters 
from 5 to 3 and as a side effect formulated a starting point for analytical 
calculations. 
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B Typical pEirameter values 

The following table contains the typical parameter values used in the simula- 
tions of lane formation in Sec. 3.2. 



Description 


Symbol 


Value 


steps to happy transition 




4 


steps to unhappy transition 




3 


floor field diffusion 


D 


0.01875 


fioor field decay 


5 


0.005 


first active floor parameter 


bi 


0.15 


second active floor parameter 


b2 


0.15 


first passive floor parameter 


9i 


0.23 


second passive floor parameter 


92 


0.10 



The first two lines show the number of consecutive allowed moves which a 
particle in unhappy mode needs to become happy again and the number of 
consecutive forbidden moves for the inverse transition, respectively. 

The following two lines give the parameters for the modification of the floor 
field as shown in (7). 

The active floor parameters determine the influence of the floor field on the 
particles (8), whereas the passive floor parameters describe the action of the 
particles on the floor field (9). 

The first six of these values can be found in the lower half of the configuration 
panel in Figure 6 (the diffusion constant is scaled by |). 
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